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Spin manipulation using electric currents is one of the most promising direc- 
tions in the field of spintronics. We used neutron scattering to observe the 
influence of an electric current on the magnetic structure in a bulk material. 
In the skyrmion lattice of MnSi, where the spins form a lattice of magnetic 
vortices similar to the vortex lattice in type II superconductors, we observe 
the rotation of the diffraction pattern in response to currents which are over 
five orders of magnitude smaller than those typically applied in experimental 
studies on current-driven magnetization dynamics in nanostructures. We at- 
tribute our observations to an extremely efficient coupUng of inhomogeneous 
spin currents to topologically stable knots in spin structures. 
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The discovery of the effect of giant magnetoresistance, now used commercially in hard disk 
drive industry, is widely recognized as the starting point of the field of spintronics. It represents 
the first example of electric currents controlled efficiently by spin structures. The comple- 
mentary process of so-called spin transfer torques, where magnetic structures and textures are 
manipulated by electric currents (i, 2), appears to be even more promising. For instance, strong 
current pulses allow to move ferromagnetic domain walls (3, 4), switch magnetic domains in 
multilayer devices (5, 6), induce microwave oscillations in nanomagnets (7) and switch ferro- 
magnetic semiconductor structures (8). However, the typical current densities required to create 
observable spin transfer torques in present day studies exceed lO^^Am"^. Because this implies 
extreme ohmic heating it was generally believed that spin torque effects can be studied exclu- 
sively in nanostructures. We report the observation of spin transfer torques in a bulk material, 
the skyrmion lattice phase of MnSi. The spin transfer torques appear when the current density 
exceeds an ultra-low threshold of ~ 10^ Am~^ - five orders of magnitude smaller than those 
used typically in experimental studies on current-driven magnetization dynamics in ferromag- 
netic metals and semiconductors. 

The skyrmion lattice in chiral magnets, like MnSi and related B20 compounds, was only 
recently discovered in neutron scattering studies (9, 10, 11, 12) and confirmed to exist in Lorentz 
force microscopy for Fei-^Co^Si {x — 0.5) (13). It represents a new form of magnetic order that 
may be viewed as a crystallization of topologically stable knots of the spin structure that shares 
remarkable similarities with the mixed state in type II superconductors. For zero magnetic field 
(Fig. 1 A) helimagnetic order appears in MnSi below Tc = 29.5 K. In a small magnetic field the 
skyrmion lattice stabilizes in a pocket below Tc, also known as the A phase. The spin structure 
of the skyrmion lattice in MnSi consists of a hexagonal lattice of magnetic vortex lines oriented 
parallel to the magnetic field B (inset to Fig. 1 A). 

Skyrmion lattices in chiral magnets are attractive for studies of spin torque effects, because 
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they are coupled very weakly to the atomic crystal structure (9) and may be expected to pin 
very weakly to disorder. In addition, electric currents couple very efficiently to skyrmions 
as follows. When the conduction electrons in a metal move across a magnetic texture, their 
spin follows the local magnetization adiabatically. Spins which change their orientation pick 
up a quantum mechanical phase, the Berry phase, that may be viewed as an Aharonov-Bohm 
phase arising from a fictitious effective field {14, 15, 16) B^g — ^eijkM{djM x dkM), where 
M = M/|M| is the direction of the local magnetization and $o — h/eis the flux quantum for 
a single electron. In the skyrmion lattice has a topologically quantized average strength of 
— $0 per area of the magnetic unit cell (for MnSi Bes « 2.5 T (11)). Bes induces an effective 
Lorentz force which gives rise to an additional "topological" contribution to the Hall effect 
proportional to the product of B^f^ and the local polarization of the conduction electrons as 
observed experimentally (11, 17). Correspondingly, because the electrons are deflected, a force 
is exerted on the magnetic structure so that there is an efficient 'gyromagnetic coupling' (18) of 
the current to the skyrmion lattice (19). 

From an alternative point of view, the skyrmion lattice may be viewed as an array of cir- 
culating dissipationless spin currents, because the skyrmions are characterized by gradients in 
the spin-orientation related to their quantized winding number. This is analogous to supercon- 
ductors, where dissipationless charge currents flow around quantized vortices due to gradients 
of the phase. When an extra spin current is induced by driving an electric current through the 
magnetic metal, the spin currents on one side of the skyrmion are enhanced while they are re- 
duced on the other side. As for a spinning tennis ball, this velocity difference gives rise to a 
Magnus force acting on the skyrmions. Note, however, that spin (due to spin-orbit coupling) 
is in contrast to charge not conserved and therefore this intuitive picture is incomplete. Most 
importantly, also further dissipative forces arise (20) which drag the skyrmions parallel to the 
current. 
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In Fig. 1 B this magnetic Magnus force, which is perpendicular to current and field direction, 
is sketched together with the additional drag forces. The Magnus and drag forces may lead 
to a translational motion of the skyrmion lattice. However, for the current densities used in 
our experiment the drift velocity of the electrons and therefore also the drift velocity of the 
skyrmions is not very large and thus very difficult to detect in a scattering experiment. 

In contrast to a translational motion, a rotation is much easier to measure in neutron scat- 
tering. Thus, we performed our experiment in the presence of a small temperature gradient 
parallel to the current, causing the magnetization and therefore the spin currents to vary in mag- 
nitude across a domain of the skyrmion lattice. In turn, the strength of the Magnus force varies 
across the skyrmion lattice (Fig. 1 B), inducing a net torque. As estimated below, the torques 
are sufficiently strong to induce rotations which can be measured directly by neutron scattering. 

For our measurements an electric current was applied along bar-shaped single crystals, 
where the direction of the current was always perpendicular to the magnetic field and there- 
fore to the skyrmion lines. In the following the neutron beam was always coUinear to the 
magnetic field (27). The six-fold diffraction pattern of the skyrmion lattice at zero current, 
j — 0, (Fig. 2 A) can be compared to the same scattering pattern at a current density, j — 
2.22 • 10^ A m~^, first in a set-up minimizing any thermal gradients along the current direction 
(Fig. 2B). The current was applied along the vertical [110] direction, whereas the field and the 
neutron beam were coUinear to the line of sight and along [110] (the horizontal direction is along 
[001]). Under current the peaks of the diffraction spots remain in the same location and broaden 
azimuthally. 

We next generated a small temperature gradient along the direction of the current as ex- 
plained in Ref. (21). Fig. 2C shows the diffraction pattern of the skyrmion lattice for this 
set-up under an electrical current density, j = 2.22 • 10® Am~^, which shows a pronounced 
counter-clockwise rotation as compared to Figs. 2 A and B (note that arrows show the technical 
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current direction and the line of sight is opposite to the direction of the neutron beam consis- 
tent with convention). When the current direction is reversed, the rotation changes sign and the 
diffraction pattern turns clockwise (Fig. 2D). 

There are several unusual aspects of this rotation. First, the entire scattering pattern rotates 
with respect to its center, i.e., all spots move by the same angle even though the electric current 
has a distinct direction. Second, when reversing either the direction of the current or the direc- 
tion of the applied field the sense of the rotation changes sign. This is illustrated in Figs. 2 E 
and F, which show the difference of intensity under current reversal and simultaneous reversal 
of current and field, respectively. For the latter case the difference of intensities vanishes. 

To confirm that the small temperature gradient along the current direction causes the rota- 
tion of the scattering pattern, we reversed the direction of the thermal gradient. As illustrated in 
Fig. 2 G and H this reverses the sense of rotation with respect to the current and field direction 
appUed in Fig. 2 C and D. Thus the differences of intensity under field reversal, shown in Fig. 2 1, 
are reversed as compared with Fig. 2 E (red and blue spots have changed location). When revers- 
ing both current and field direction, the difference of the patterns. Fig. 2 J, vanishes as before. 
Finally, when applying the current along a different crystallographic direction (we tested (111) 
and an arbitrarily cut sample) the same antisymmetric rotations of the diffraction pattern as a 
function of magnetic field, electric current and temperature gradient are observed {21). 

The detailed rotation as a function of the applied current was determined in systematic 
measurements for various temperatures and samples. Here the temperature measured at a spe- 
cific spot at the surface of the sample {21) was kept constant and temperature gradients had 
always the same direction. The resulting current dependence of the azimuthal rotation angle 
A$, shown for three temperatures as a function of the applied current density in Fig. 3, exhibits 
a well-defined threshold of order jc ~ lO^A above which the rotation begins. For \ j\ > jc 
the entire scattering pattern rotates and the rotation angle A$ increases steeply with increasing 
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Potential parasitic effects cannot explain the observed rotations. First, detailed studies show 
that there are no changes of orientation of the skyrmion lattice as a function of temperature 
at j = (10, 12, 13). Second, the temperature difference AT between the sample surface 
and the sample support shows a smooth quadratic increase with current density independent 
of the direction of the current (Fig. 3 B) in contrast to Fig. 3 A. Third, for the current densities 
appUed in our study the Oersted field increases from zero (at the center of the sample) to a 
value of roughly 1 mT at the surface of the sample - much smaller than the applied magnetic 
field of 175 mT and therefore negligible. Finally, for a current parallel to the skyrmion Unes or 
the pristine helimagnetic state neither a rotation nor a broadening are observed. Interestingly, 
recent numerical simulations (22, 23) suggest that much larger current densities of the order 
10^^ A m~^ may change the orientation of helical magnetic structures. 

To explain our experiments the interplay of three tiny forces has to be considered: (i) spin 
transfer torques, i.e., current induced forces, (ii) pinning forces and (iii) anisotropy terms. These 
determine the origin of the rotation, the presence of a threshold and the angle A$ of the rotated 
diffraction pattern, respectively. The spin transfer torques can, e.g., be described by a Landau- 
Lifshitz-Gilbert (LLG) equation or variants of Landau-Lifshitz Bloch equations (24) which in- 
clude both reactive and dissipative components representing the Magnus and drag forces men- 
tioned above, respectively. Both are expected to be of the same order of magnitude (25). While 
the strength of the dissipative forces (related to the parameter ^ in the LLG equations) is not 
known, it is possible to estimate the strength of the reactive forces quantitatively. 

The size of the Magnus force is given by the product of spin current and the fictitious 
effective magnetic field, /m ~ ejgBes, and may be estimated as 
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where a fa 4.58 A is the lattice constant of MnSi, Tc 29.5 K the ordering temperature, and the 
local temperature dependent polarization is defined as the ratio of the spin and charge current 
densities times the elementary charge, p{T) — ejs/j. For the skyrmion phase we estimate 
p{T) 0.1 (11). The resulting forces at the current densities of 10^ Am~^ studied are much 
larger than, e.g., gravitational forces on the sample, but small when expressed in terms of the 
microscopic units, ksTc/a'^ (cf. Eq. 1) raising the question why the critical currents jc are so 
small. 

For j < jc the current induced forces are balanced by pinning forces caused by disorder 
and the underlying regular atomic crystal lattice. The latter may be neglected because of the 
small spin-orbit interaction (26). The discussion presented in (21) suggests, that even a very 
strong defect, which locally destroys the magnetization completely, will result in a very small 
pinning force, less than a few 10~^ ksTc/a per impurity, mainly because the magnetization of 
the skyrmions varies very smoothly. Therefore the observed critical current density jc together 
with our estimate, Eq. 1, is consistent with strong pinning defects with a density below 1 ppm. 
In fact, even though the real density of defects may be higher, their influence may be strongly 
reduced as the system is in the 'collective pinning' regime known, e.g., from vortices in super- 
conductors (27). Here pinning forces of random orientation average out to a large extent due to 
the rigidity of the skyrmion lattice. 

The size of the rotation of the skyrmion lattice for j > jc reflects the balance of the torques 
tm and tl due to inhomogeneous Magnus forces and the atomic lattice, respectively. We start 
by noting, that by symmetry the orientation of a perfect skyrmion lattice (described by a third 
rank tensor) cannot couple linearly to the current j, because of the sixfold rotational symmetry 
of the lattice. However, a small temperature gradient breaks this symmetry and generates sizable 
variations of the amplitude of the magnetization and associated polarization p(T) of the electric 
currents, because the skyrmion phase is only stable close to Tc (cf. Fig. 1). As a result, the 
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Magnus force, Eq. 1, which is proportional to the spin currents and therefore to p{T), will be 
considerably larger at the 'cold' side of the skyrmion lattice than at its 'hot' side. This gives 
rise to a net torque per volume Tm ~ /r x fM(r)ci^r/\^. Using Eq. 1 with a spatially varying 
p{T) and ignoring again dissipative forces for an order of magnitude estimate, we find that the 
rotational torque per volume, tm, in the direction of the fictitious field B^s for a skyrmion 
lattice domain of size R is given by 

where we assume j ?s jc and Vp 0.1/10 mm 

Note that the effect is proportional to the gradient of the temperature parallel to the cur- 
rent as Vp OT^-^ Thus Tm changes sign when either current, magnetic field or 
temperature gradient are reversed. The sign of tm and all other forces obtained from this anal- 
ysis (Fig. IB) is consistent with all our experiments taking into accout that charge carriers are 
hole-like as measured experimentally (11,17). This explains our main experimental results. Ac- 
cording to this analysis the rotational torques arise from temperature-gradient induced gradients 
in the spin-current. Interestingly, it is more difficult to generate an analogous effect for vortices 
in superconductors, because charge, as opposed to the spin in skyrmion lattices, is exactly con- 
served (in (28) a rotation of a superconducting vortex lattice has been induced with a bespoke 
current distribution). 

For an estimate of the factor (i?/l mm)^ in Eq. 2 a lower limit, R > 1/xm may be inferred 
from the resolution-limited rocking width of the magnetic Bragg peaks in the skyrmion phase 
when avoiding demagnetization fields (21). Yet, even a small torque tm may lead to large 
rotation angles, because the balancing torque, r^, which orients the skyrmion lattice relative to 
the atomic lattice, is tiny. Only anisotropy terms arising in high power of the spin orbit coupling 
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Xso contribute to the torque per volume tl which we estimate as (21) 

Ti~-10-2A^5O^sin(6$). (3) 

For small rotation angles the torque tl grows linearly in the rotation angle However, in 
contrast to the torques arising from the inhomogeneous Magnus force, tl is independent of the 
size R of the domains. The rotation angle $ is finally determined by the balance of tm and 
Tl. Because of the small prefactor in Eq. 3, 10~^A|o, the large rotation angles observed in our 
experiments can be explained even for moderately large domains. 

It is likely that for j > jc not only a rotation by an angle sets in but also a linear motion of 
the magnetic structure, because any rotation of a sizable magnetic domain requires a depinning 
from defects. For moving domains, spin currents in a frame of reference that is comoving with 
the domain enter in all formulas given above. Therefore the size of the rotation in the end also 
depends sensitively on the frictional forces which break Galilean invariance (25). 

Our observations identify chiral magnets and systems with nontrivial topological properties 
as ideal systems to advance the general understanding of the effects of spin transfer torques. For 
instance, spin transfer torques may even be used to manipulate individual skyrmions, recently 
observed directly in thin samples (13). In fact, even complex magnetic structures at surfaces 
and interfaces may be expected to exhibit the spin torque effects we report here (29, 30). 
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Figure 1 : (A) Magnetic phase diagram of MnSi. The inset schematically shows the spin struc- 
ture of the skyrmion lattice in a plane perpendicular to the applied field. (B) Schematic depic- 
tion of the spin transfer torque effects on the skyrmion lattice. A temperature gradient induces 
inhomogeneous Magnus and drag forces and therefore a rotational torque. 
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Figure 3: (A) Change of the azimuthal angle of rotation of the scattering pattern as a function 
of current density for three different temperatures. Data were recorded in a magnetic field of 
0.175 T. Above a clear threshold of 10^ Am~^ an increasingly strong rotation is observed, where 
the sign of the rotation depends on the direction of the current. (B) Temperature difference 
between the surface of the sample and the sample holder as a function of current density, where 
the temperature at the surface of the sample was kept constant for each of the three values given 
in Fig.3A. 
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1 Experimental Methods 

1.1 Neutron Scattering 

Our neutron scattering measurements were performed at the diffractometer MIRA at FRM II at 
the Technische Universitat Miinchen (SI). Data were recorded for an incident neutron wave- 
length A — 9.6 A with a 5% FWHM wavelength spread. A delayline ^He area detector of 
200 X 200 mm^ was used with a position resolution of order 2x2 mm^. 

The incident neutron beam was colUmated over a distance of 1.5 m using apertures 4 x 
4 mm^ and 1.5 (width) x 3 (height) mm^ after the monochromator and before the sample, 
respectively. The size of the sample was 1.5 (width) x 8 (height) mm^, where we used a Cd 
aperture of 3 x 3 mm^ directly at the sample in most experiments. The distance from the sample 
to the detector was between 0.8 and 1.3 m. 

1.2 Cryogenic environment and magnetic field 

Two different orientations between magnetic field, current direction and neutron beam were 
used. In set-up 1 the neutron beam was parallel to the applied magnetic field (red arrow) that 
stabilized the skyrmion lattice. The applied electric current (green arrow) was perpendicular 
to the neutron beam and the magnetic field. In set-up 1 the current was flowing perpendicular 
to the skyrmion lines. In set-up 2 the incident neutron beam was perpendicular to the applied 
magnetic field and the current. Therefore, in set-up 2 the current was flowing parallel to the 
skyrmion lines. 

The magnetic field was generated with bespoke water-cooled Cu solenoids in a Helmholtz 
configuration (S2). The magnetic field profile was carefully characterized with a Hall probe and 
found to be uniform better than 1% over the sample volume. In the A-phase data were identical 
after zero-field cooling and field cooling. 
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Figure S 1 : (A) Set-up 1 used for all data shown in the main text. The incident neutron beam 
was parallel to the applied magnetic field (red arrow) that stabilized the skyrmion lattice, i.e. 
the neutron beam was parallel to the skyrmion lines. The applied electric current (green arrow) 
was perpendicular to the neutron beam and the magnetic field and hence skyrmion lines. (B) 
Set-up 2, in which the incident neutron beam was perpendicular to the applied magnetic field 
and the current. In this set up the current was flowing parallel to the skyrmion lines. 
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1.3 Sample holder 

For our measurements we used a bespoke sample stick that permitted application of currents up 
to several 10 A (55). The sample stick and sample holder were cooled with a cryogen-free pulse 
tube cooler, where He exchange gas with a pressure of several hundred mbar and spring-loaded 
mechanical contacts were used to ensure excellent thermal coupling. 

Two different sample holders were used as shown in Figs. S2 and S3. The salient features 
reported in our paper were observed with both sample holders A and B. Sample holder A con- 
sisted of a large Cu heat sink, to which the sample and current leads were attached. Most of our 
data were measured with sample holder A, in particular all the data shown in the main text. 

To insulate the sample electrically from the heat sink we used cigarette paper and a thin 
layer of GE varnish. In the low temperature section of the sample stick and sample holder the 
current leads were made of electrolytic (high-purity) Cu, soft-soldered directly to the samples 
with the help of special solder flux. 

We generated a small temperature gradient along the direction of current flow with a small 
wedge of GE varnish between the sample and the sample holder. This way the thermal coupling 
to the sample holder changed along the sample, causing a temperature gradient along the direc- 
tion of the electric current in the presence of the resistive heating by the sample. The tiny tilt of 
the sample caused by the GE varnish has no observable effect as we checked explicitly. 

In both sample holders a Pt sensor was attached to the sample surface (the sensor did not 
generate any noticeable background signal in neutron scattering). Additional temperature sen- 
sors, e.g., as shown in Fig. S2, were used to monitor the temperature difference between sample 
surface and heat sink (main body of the sample holder). 

Typical temperature differences between the surface of the sample and the heat sink for the 
case of sample holder A are shown in Fig. 3 (B) of the main text. These data were obtained for 
fixed sample temperature while the temperature of the heat sink was adjusted. With increasing 
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Figure S2: Sample holder A used in our experiments. The sample and current leads were 
connected to a large Cu heat sink. Data shown in the main text were recorded with this sample 
holder. The sample holder was cooled with a cryogen-free pulse tube cooler, where He exchange 
gas with a pressure of several hundred mbar and spring-loaded contacts were used for thermal 
coupling. 

sample holder B sample thermometer thermal anchor 




current leads 



Figure S3: Sample holder B used in our experiments. The sample holder was cooled with a 
croygen-free pulse tube cooler, where He exchange gas with a pressure of several hundred mbar 
and spring-loaded contacts were used for thermal coupling. 
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current density the temperature of the sample increases as the square of the current density, i.e., 
the increase is very smooth. This smooth increase contrasts the abrupt onset of the rotation 
above the critical current density jc as discussed in the main text. 

Fig. S4 (A) shows sample (i) as firmly attached to the heat sink. Data shown in Fig. 2 (B) of 
the main text were recorded with this configuration. Fig. S4 (B) shows sample (i) as attached 
to the heat sink with a small wedge of GE varnish. To obtain a well-defined angle a small 
support (white) was placed underneath the sample. The cold spot is at the right hand side, which 
corresponds to the top of Fig. 2 of the main text (thermal gradient pointing down). Data shown 
in Fig. 2 (C) through (F) of the main text were recorded with this configuration. Fig. S4 (C) 
shows sample (i) as attached to the heat sink, where the direction of the wedge of GE varnish is 
reversed. The cold spot is at the left hand side, which corresponds to the bottom of Fig. 2 of the 
main text (thermal gradient pointing up). Data shown in Fig. 2 (G) through (J) of the main text 
were recorded with this configuration. 

1.4 Samples 

We studied six different samples as summarized in Table SI. The samples were either prepared 
by Bridgman growth or optical float-zoning. All samples studied were single crystals with high 
residual resistance ratios (RRR) around 100 and residual resistivity around po ~ 2 fiilcm. This 
corresponds to mean free paths up to 1000 A (S4). The magnetic and thermodynamic properties 
of samples cut from the same ingots agreed very well with the literature, i.e., all samples were 
of high quality. 

For the experiments the samples were oriented with Laue x-ray diffraction. Note that data 
for the same samples has also been reported in the other publications, e.g., sample (ii) corre- 
sponds to "sample 2" in Ref. (55). We note in particular, that sample (v) exclusively served to 
study the role of demagnetizing fields as described below. All samples showed the same salient 
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sample thermometer 
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Figure S4: (A) Sample (i) as firmly attached to the heat sink. Data shown in Fig. 2 (A) and 2 (B) 
of the main text were recorded with this configuration. (B) Sample (i) as attached to the heat 
sink with a small wedge of GE varnish. To obtain a well-defined angle a small support (white) 
was placed underneath the sample. The cold spot is at the right hand side, which corresponds 
to the top of Fig. 2 of the main text. Data shown in Fig. 2 (C) through (F) of the main text were 
recorded with this configuration. (C) Sample (i) as attached to the heat sink, where the wedge of 
GE varnish is inverted. The cold spot is at the left hand side, which corresponds to the bottom 
of Fig. 2 of the main text. Data shown in Fig. 2 (G) through (J) of the main text were recorded 
with this configuration. 
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Figure S5: (A) Scattering pattern of sample (vi) without current and hence temperature gradient. 
(B) Scattering pattern of sample (vi) under current and small temperature gradient. (C) same as 
panel B with current direction reversed. 

properties under applied currents. This includes in particular sample (iii) which had a different 
aspect ratio of the cross-section and sample (vi), in which the current was applied along (111). 
Typical data recorded in sample (vi) are shown in Fig. S5. In other words, the results reported 
in the main text do not depend on the specific crystallographic orientation and sample shape. 
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Table SI: MnSi single crystals investigated in this study. 'Length' refers to the size perpen- 
dicular (vertical) to the beam direction and represents the direction of current flow; 'thickness' 
refers to the size in the direction of the neutron beam; 'width' refers to the size perpendicular 
(horizontal) to the beam direction. 
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1.5 Data for currents applied parallel to the Skyrmion lattice 

Typical data of the skyrmion lattice for the applied magnetic field perpendicular to the incident 
neutron beam (set-up II in Fig. SI) are shown in Fig. S6. Data were recorded for sample (ii), 
which showed a slightly lower critical current than sample (i). Data were recorded with sample 
holder B, where Fig. S6 (A) shows the scattering pattern for j = 0. 

Due to a small temperature gradient that exists in sample holder B at j = there are weak 
diffraction spots of the conical phase (spots at bottom and top) even at j — 0. The horizontal 
spots labelled '3' and '6' correspond to spots '3' and '6' of the six-fold patterns shown in Fig. 2 
of the main text (data in Fig. 2 were recorded in sample (i)). The diffraction pattern under an 
appUed electric current density of j — 1.24 • 10^ A/m^ parallel to the skyrmion lines is shown in 
Fig. S6 (B). Neither a rotation nor broadening may be seen. Note that the same sample (number 
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Figure S6: (A) Diffraction pattern in the skyrmion lattice phase for magnetic field perpendicular 
to the incident neutron beam. Data were recorded in sample (ii). The spots labelled as '3' and 
'6' correspond to the same spots shown in Fig. 2 of the main text for magnetic field parallel to 
the neutron beam. (B) Diffraction pattern in the skyrmion lattice for electric current parallel 
to the skyrmion lines. Note that the spots of the skyrmion lattice do not show any rotation or 
broadening whatsoever. For the same current density applied perpendicular to the skyrmion 
lines strong rotations are seen. 

(ii)) as studied on the same sample holder shows strong rotations and broadening when the same 
current is applied perpendicular to the skyrmion lines just like the features reported in the main 
text (cf. Fig. 2). 

1.6 Data Analysis 

Neutron data were analyzed with the software GRASP, v4.26, developed at the ILL. For each 
data set recorded at fixed current and temperature we determined at first accurately the center 
of the scattering pattern. We found slight systematic variations of the center of the scattering 
pattern as a function of current density. All azimuthal positions and widths of the intensity 
maxima reported in our manuscript were determined with respect to the center of the pattern. 
Data were binned as a function of azimuthal angle in intervals of 1 degree in a ring with respect 
to the center of the scattering pattern encompassing the scattering maxima. The azimuthal 
position and width was determined by fitting a simple Gaussian function. 
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Figure S7: Schematic view as seen from above of the variation of the magnetic field lines across 
the sample. Due to demagnetizing effects at the fringes of the sample the internal field neither 
coincides with the tilting angle nor with the incident neutron beam. 

1.7 Role of Demagnetizing Fields 

The mosaic spread of the skyrmion lattice was inferred from rocking scans. In some of the 
experiments we could only perform rocking scans with respect to a vertical axis. As a result 
the spots at the top and bottom may be much weaker because they do not satisfy the scattering 
condition. 

In this section we summarize a recent study, where we compared the apparent mosaic spread 
derived from the rocking width of the magnetic satellites in the skyrmion lattice for different 
sample cross-sections. A full account of this work will be reported elsewhere (5(5). In our 
study we find that the hexagonal magnetic scattering intensity in the A phase aligns strictly 
perpendicular to the applied magnetic field. This implies that the magnetic scattering intensity 
depends sensitively on variations of the magnetic field directions across the sample volume, 
notably the effects of demagnetizing fields. 
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An easy way to notice the influence of the demagnetizing fields was in the context of the 
orientation of the hexagonal scattering intensity with respect to the orientation of the sample as 
described in the supplementary online information of Ref. (S5). When the appUed field was not 
perfectly perpendicular to the disc, we observed a deflection of the maxima consistent with the 
demagnetizing fields as illustrated in Fig. S7. 

In some of our studies rocking scans could only be carried out for the vertical axis, i.e., along 
the long side of our sample which was parallel to a (110) direction. It was oriented vertical and 
perpendicular to the incident neutron beam and magnetic field as shown in Fig. SI (A). Thus in 
the A-phase two spots of the hexagonal scattering intensity coincided always with the vertical 
(110) axis and thus the axis of the rocking scans. 

For a reasonably well coUimated neutron beam the observation of equal intensity of the 
vertical spots implies, that neither of these vertical spots satisfies the Bragg condition. The in- 
tensity originates only from the tails of the Bragg spot. Therefore, when summing over rocking 
scans with respect to the vertical axis vertical spots remain weak. 

Shown in Fig. S8 are the rocking scans for sample (ii) studied in the spin torque measure- 
ments. In the helical state the half- width of the rocking scans corresponded to a magnetic mo- 
saicity 1]^ ^ 3.5° consistent with previous work and long range order (57). In the A-phase the 
half- width of the rocking scans corresponded to an apparent magnetic mosaicity rjm ~ 1.75°. 

In order to identify the origin of the peculiar angular dependence observed in the skyrmion 
lattice in the spin torque sample we measured rocking scans of the helical order and the skyrmion 
lattice in a disc that was only a few tenths of a mm thick with a diameter of roughly 15 mm. 
The disc was positioned perpendicular to the magnetic field and incident neutron beam. Using 
a small source aperture only a tiny cross-section of a few mm^ in the center of the disc was 
studied, thereby avoiding any anomalous internal field distributions at the fringes of the sample. 

Typical rocking scans observed in this configuration are shown in Fig. S9. For convenience 
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rocking angle {") 

Figure S8: (A) Rocking scans in the helical state of MnSi. A Gaussian distribution with a width 
at half-maximum corresponding to a magnetic mosaicity r/^ ~ 3.5° is observed consistent with 
previous studies and the literature (57). (B) Rocking scans in the A-phase of MnSi. The ap- 
parent width is somewhat reduced, however, the functional dependence is completely different. 
As shown in Fig. S9 this distribution is most likely due to demagnetizing effects. See also the 
onhne supplementary information of Ref. (55). 
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Figure S9: Typical intensity variation in rocking scans in the helical state and the skyrmion 
lattice of MnSi for a sample (v). Only the central section of the sample was illuminated by the 
neutrons, thereby avoiding any fringe effects at the rim of the sample. This way the internal 
field distribution was extremely narrow. The rocking scans show a width Aco ~ 3.5° at half 
maximum in excellent agreement with all previous studies and the literature. In contrast, in the 
skyrmion lattice the rocking width is resolution limited Ao; ~ 0.4°. This suggests strongly that 
the skyrmion lattice is very well ordered and rather rigid as assumed in the discussion in the 
main text. 
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the intensity is normalized. In the helical state the rocking width is Aoo ^ 3.5°, in excellent 
agreement with values reported in the literature (57). This value also agrees very well with the 
rocking width observed in the spin torque sample. Thus the rocking width appears to reflect the 
intrinsic magnetic mosaicity of the helimagnetic state in MnSi. 

In contrast, the rocking width in the skyrmion lattice is now resolution limited with Aw ~ 
0.4° - almost an order of magnitude smaller. Astonishingly, the corresponding coherence length 
of the magnetic state is of the order of several ten micro-meters ! In turn this observation strongly 
suggests, that the rocking width observed in the spin torque sample shown in Fig. S8(B) results 
from the distribution of demagnetizing fields in the sample. In other words, the skyrmion lines 
closely follow the internal field distribution and the skyrmion lattice exhibits well developed 
long-range order. This is an important assumption for the analysis of our data in the main text. 
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2 Theoretical analysis 

The magnetic properties of MnSi are well described by a Ginzburg-Landau theory for the local 

magnetization M(r), 

F[M] = j d^r (roM^ + J{VM.f + 2L) M ■ (V x M) + UM.'^ - B ■ m) , (SI) 

where B is the external magnetic field and tq, D, J,U are parameters (U, J > 0) (S8, S9). As 
shown in Ref. (55), this theory describes the formation and phase diagram of the skyrmion 
lattice if thermal fluctuations around the mean field solution, i.e. around local minima of (SI), 
are taken into account. The energy scales and the hierarchy of forces in the system are governed 
by the strength of spin-orbit coupling, \so, which is very small in MnSi. The most important 
component attributed to spin-orbit coupling is thereby the Dzyaloshinsky Moriya interaction, 
D oc \so, which determines the size and distance of the skyrmions given hy J/D oc l/Xso 
(for details see Ref. (55)). 

2.1 Orientation of the Skyrmion lattice 

The Ginzburg-Landau potential in lowest order in Xso given in Eq. (S 1) is rotationally invariant 
around the axis defined by the magnetic field B. The orientation of the skyrmion lattice within 
the plane perpendicular to B in the absence of a current, j — 0, is determined by the anisotropy 
of the atomic crystal, i.e., anisotropy terms of higher order in Xso that are not yet included in 
Eq. (SI). Such anisotropics favor certain orientation angles $ of the magnetic texture. Due to 
the sixfold symmetry of the skyrmion lattice, a ^-potential proportional to at least — cos(6$) 
is needed in order to orient it. Such anisotropy terms are, however, suppressed by the tiny 
factor A|o with Xso ~ 10"^. An example of such a term is Jid^Mf + {d^Mf + {d^Mf 
which is proportional to Xgo because the distance between the skyrmions, ~ J/ D, is linear in 
I /Xso- In contrast, the energy density, e^, of the skyrmion lattice, which according to specific 
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heat measurements (SIO) is of the order of AF/V ~ IO^^/cb^'c/ o^, is only proportional to 
A|o (55). The torque per volume resulting from such a potential can therefore be estimated as 
tl ~ A|o9$ cos(6<l>) AF/\/, which reproduces Eq. (3) in the main text. 

2.2 Disorder pinning of the Skyrmion lattice 

In our spin torque experiment an unexpectedly small threshold current jc ~ 10~^ Am^^ was 
observed (Fig. 3(A) of the main text) suggesting an inefficient pinning of the Skyrmion lattice 
by disorder. It is instructive to compare this observation with the seemingly similar situation 
encountered in type II superconductors where the pinning of vortex lines has been studied ex- 
tensively (Sll). In superconductors the order parameter varies on a length scale typically much 
smaller than the distance between the vortices and, in particular, vanishes in the vortex cores. A 
pinning force F due to an impurity can be estimated by taking the energy difference AE of a 
defect in the core center and a defect outside of the superconducting vortex core divided by the 
typical radius, Tc, of the core, F ~ AE/tc- 

In contrast, pinning may be expected to differ significantly in the case of a Skyrmion lattice 
because the magnetization is smooth and does not vanish anywhere (85). Besides the low defect 
concentration of our samples, evident from the large charge carrier mean free paths around 
1000 A, there are several aspects that have to be taken into account when estimating the pinning 
force arising from a single, localized defect: (i) The free energy density of the skyrmion lattice 
as estimated from the measured specific heat {SIO), AF/V ~ lQ~'^kBTc/a^, is tiny, (ii) As 
already mentioned above, pinning forces arising from the coupling to the magnitude of the 
magnetization are small because it varies by less than ±20% (55) within the magnetic unit cell, 
(iii) The coupUng of disorder to the direction of the spins originates from spin-orbit interactions. 
It is hence also very small, (iv) The forces, i.e., the energy changes per length, are small 
because the distance between the skyrmions is a factor 40 larger than a. Combining all these 
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factors, we estimate the pinning force from a single strong defect (e.g. by replacing a Mn 
atom by a nonmagnetic impurity) to be less than a few lO^^/c^Tc/a. Here we assumed that 
the magnetization is completely destroyed in a fraction of the unit cell (there are 4 Mn and 4 
Si atoms per unit cell) and that the difference of the free energy for a defect in the center of a 
skyrmion and, e.g., between skyrmions is only a fraction of /\F/V due to the small variations 
of the amplitude of the magnetization. The most important factor is, however, the large distance 
between the skyrmions. 

As discussed in the main text, the total pinning force per volume can not be obtained by 
multiplying the pinning force of a single defect by the defect density, because the pinning 
forces from several defects within a domain can partially average out if the skyrmion lattice 
is sufficiently rigid. The analogous collective pinning regime has been studied extensively for 
superconducting vortices {Sll). For a full quantitative theory of pinning, one would therefore 
need information on (i) the nature of single pinning centers, (ii) the defect concentration and 
(iii) the rigidity of the skyrmion lattice. While the latter can be calculated or measured with 
neutron scattering and while some information of defect concentration can in principle be in- 
ferred from the residual resistivity or de-Haas-van- Alphen measurements, the nature of defects 
important for either magnetic pinning or scattering of electrons is presently not clear. Neverthe- 
less, we conclude that besides the small defect concentration the smoothness of the skyrmion 
spin structure is probably the most important factor to explain the small pinning forces. 
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